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Abstract: We revisit the issue of probing R- violating couplings of supersymmetric 
theories at hadronic colliders, particularly at the LHC. Concentrating on dimuon 
production, an evaluation of the optimal sensitivity to the i?-violating coupling is 
performed through a maximum likelihood analysis. The measurement uncertainties 
are evaluated through a study of fully generated events processed through a fast 
simulation of the ATLAS detector. It is found that a host of .R-violating couplings 
can be measured to a statistical accuracy of better than 10%, over a significant part 
of the mj - A parameter space still allowed by low energy measurements. Since the 
bounds thus obtained do not simply scale as the squark mass, one can do significantly 
better at the LHC than at the Tevatron. The same analysis can also be extended 
to assess the reach of the LHC to effects due to any non-SM structure of the four- 
fermion amplitude, caused by exchanges of new particles with different spins such as 
leptoquarks and gravitons that are suggested by various theoretical ideas. 
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1. Introduction 

Within the Standard Model (SM), electroweak gauge invariance ensures the conser- 
vation of both lepton number and baryon number, at least in the perturbative con- 
text. However, this is not so within the Minimal Supersymmetric Standard Model 
(MSSM). The most general superpotential respecting the gauge symmetries of the 
SM contains bilinear and trilinear terms which do not conserve either the baryon 
number (B) or the lepton number (L). Clearly, the simultaneous presence of both 
lepton- and baryon-number-violating operators could lead to very rapid proton de- 
cay, especially for TeV scale sparticle masses. The existence of all such terms can be 
forbidden by postulating a discrete symmetry [1,2], called .R-parity, which implies a 
conserved quantum number R p = (_i) 3B + L + s ; where 5* stands for the spin of the 
particle. The very definition implies that all the SM particles have R p = +1, while 
all the superpartners are odd under this symmetry. Thus, apart from suppressing 
proton decay, it also guarantees the stability of the lightest supersymmetric particle 
(LSP), thereby offering a ready-made candidate for cold dark matter. 

However, while a conserved i?-parity seems desirable, it is perhaps too strong a 
requirement to be imposed. For one, the measure is an ad hoc one and there does 
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not exist an overriding theoretical motivation for imposing this symmetry, especially 
since a suppression of the proton decay rate could as well be achieved by ensuring that 
one of B and L is conserved. Indeed, it has been argued [3] that this goal is better 
served by imposing a generalized baryon parity instead. Unlike -R-parity, this latter 
(Z3) symmetry also serves to eliminate dimension-5 operators that could potentially 
have led to proton decay. Furthermore, non-zero lfi p couplings provide a means of 
generating the small neutrino masses that the neutrino oscillation experiments seem 
to call for. Similarly, a significant value for such couplings has been shown to pro- 
vide respite from the tachyonic nature of sleptons in models with anomaly-mediated 
supersymmetry breaking [4]. It is thus of both theoretical and phenomenological 
interest to consider violations of i?-parity. 

Limiting ourselves to a renormalizable superpotential, the possible i? p -violating 
terms can be parametrized as 

W D ^H 2 + E fakLiLjEZ + K jk LiQjD c k + K^D)D^ ) (1.1) 

i i,j,k 

where i,j,k are generation indices, L (Q) denote the left-handed lepton (quark) 
superfields, and E c , D c and U c are the right-handed superfields for charged leptons, 
down and up-type quarks respectively. The couplings \jk and \"j k are antisymmetric 
in the first and the last two indices, respectively. A conserved baryon number requires 
that all the \'^ k vanish identically, thereby avoiding rapid proton decay. 

As with their i? p -conserving cousins, namely the usual Yukawa couplings, these 
couplings are entirely arbitrary. Some phenomenological constraints exist, though. 
For example, the preservation of a GUT-generated B — L asymmetry necessitates 
the preservation of at least one of the individual lepton numbers over cosmological 
time scales [5]. At a more prosaic level, the failure of various collider experiments [6] 
to find any evidence of supersymmetry 1 has implied constraints in the parameter 
space. Even for superpartners too heavy to be produced directly, strong bounds 
on these couplings can be inferred from the remarkable agreement between low en- 
ergy observables and the SM predictions. These include, for example, meson decay 
widths [9,10], neutrino masses [10,11], rates for neutrinoless double beta decay [12], 
etc. The bounds generally scale with the sfermion mass and, for mj = 100 GeV, they 
range from ~ 0.02 to 0.8 [13]. 

In general, therefore, when one studies the collider signals for i? p -violating su- 
persymmetry, it is usual to consider the so-called 'weak' i? p -violation scenario where 
the production of the superparticles goes through gauge couplings and the only role 
of R p - violation is in the decay of the lightest supersymmetric particle [14]. Such 
studies are clearly insensitive to the exact size of the i? p -violating coupling as long as 

1 Although i?-parity violation has been touted as an explanation [7] for the reported excess of 
high-Q 2 events at hera [8], it is no longer clear that this anomaly persists. 
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it is large enough to make the decay length of the LSP undetectable 2 . The processes 
that are directly sensitive to the size of this i? p -violating coupling are the produc- 
tion of sparticles through them [7, 16], the decays of sparticles through them [17, 18] 
or through indirect effects of the virtual sparticle exchange by interference of the 
R p - violating amplitude with the SM one [17,19-22]. 

An inspection of the aforementioned low-energy constraints shows that they are 
strongest when the term involves only the first two generations, and are often rather 
weak when one or more of the superfields belong to the third generation. In the 
context of collider experiments, a measurement of such couplings could, in principle, 
be done in more than one way, e.g., direct production, the decays of particles such 
as the top quark or sparticles other than the LSP through coupling, as well 
as through virtual exchanges. For example, in the event of a large coupling, one 
could study the rate for the production of a single superparticle. However, such 
measurements are subject to uncertainties from the luminosity measurement in some 
cases and also from the knowledge of branching ratios. We therefore revisit the issue 
of the indirect determination of some of the A- jfc couplings at the LHC, using their 
contribution to the production of lepton pairs. 

Clearly, for such contributions to be significant, the initial state must involve 
quarks of the first generation or, in other words, at least one of j and k must be 1. 
Restricting ourselves, for the sake of definiteness, to muons in the final state, one 
sees that the least constrained among the relevant couplings [13] are A' 231 and X' 2ll . 
We look at the reach of the LHC in the two-dimensional X'-rriq plane. The current 
study goes much beyond the earlier analyses of this process [19,21] at the Tevatron 
in that we exploit the differences in the angular distribution of the leptons as well 
as the invariant mass distribution of the lepton pair. The use of this additional in- 
formation increases the sensitivity to the _R p -violating couplings. Through a detailed 
experimental analysis, we show that it should be possible to keep the systematic 
uncertainties at the level of a few per cent. We further find that this analysis would 
serve a role complementary to that played by looking for ^-violating decays of the 
squarks. 

The rest of the article is planned as follows. In the next section, we first outline 
the details of the production of pairs at the LHC along with the discriminatory 
features of the various distributions we use. We present the result of a parton- 
level calculation in this discussion. Following it up are the details of the maximum 
likelihood analysis technique used by us, and we present the results of the same at the 
parton level as well as at the level of fully generated events. We then give a discussion 
of the various systematic uncertainties in the analysis. After this, we briefly discuss 
the i?p-violating branching ratios of the third-generation squark. We show that for 
the squark mass range considered in this paper, there exists a region in parameter 

2 If any of the i? p -violating couplings is > 10~ 6 or so, then the LSP will decay within the 
detector [15]. 
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space in which the coupling could in principle be measured at the 5% level at 
LHC at full luminosity. This demonstrates the above mentioned complementarity 
very nicely. 



2. Drell— Yan in the presence of R p violation 



The Drell- Yan process at a hadronic collider has, over the years, been studied in 
great detail. Its analysis has served as an excellent theoretical laboratory, first for the 
quark parton model and later for perturbative QCD. In fact, not only the NLO, but 
even the NNLO corrections have been computed [23]. In the context of the present 
analysis, though, we shall limit ourselves to leading-order calculations. The reasons 
are manifold. For one, the NNLO calculations, within the SM, have demonstrated 
that a constant .ff -factor gives a very good description of the corrections over a very 
wide range of kinematic variables. Secondly, the radiative corrections in the presence 
of non-SM physics, such as the case under study, have not yet been calculated. It 
might be argued, though, that since we are interested in probing typically small 
$p couplings, a very precise calculation of the strong corrections to these additional 
contributions is not crucial. Instead, we might assume that the .fT-factor is essentially 
the same as in the SM. Clearly, this estimate is quite a reasonable one; it is supported 
by the recent calculations [24] of i^-factors for chargino or neutralino pair-production, 
each of which receives t-channel contributions as in the present case. This assumption 
simplifies the analysis, as the Born term is extremely simple to analyse. 

A look at the superpotential given by Eq. (|1.1[) tells us that it is only the \'^ k 
couplings that can affect Drell- Yan production of a dilepton pair. Expressed in terms 
of the component fields, the relevant part of the Lagrangian reads 



Cy 



A ijk 



-s L d k R d{ + # L db>i + (4r(^) c 4 

-el44-44ei-(4)*(ei) c <' 



+ h.c. 



(2-1) 



A non-zero A' 2jA . would then lead to an additional w-channel (d R exchange) diagram for 
the process UjUj — > and a t-channel (u J L exchange) diagram for d^dh — > 

Note that neither resonance production processes leading to the dimuon final state 
nor processes such as q^qj — > e~fi + can occur when only a single Iji p coupling is 
non-zero, an assumption [25] that we shall work with. 

The differential cross section is modified in a straightforward way and reads, in 
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Figure 1: Relative change of the dimuon cross section at the LHC as function of the 
coupling X^u- The assumed lower limit on the dilepton invariant mass is 500 GeV. The 
curves correspond to four different squark masses: 300, 500, 700 and 900 Ge V respectively. 
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where represents the charge of the fermion / and g{ R its couplings to the Z. 
Clearly, the new contribution is relevant only if q — u, d or, in other words, only 
if at least one of j and k refers to the first generation. Of the nine T^ v couplings 
that could, in principle, contribute to dimuon production, we thus need to concern 
ourselves with only five. For a given strength of the coupling, the change wrought in 
the total cross section would depend on the mass of the squark(s) involved as well 
as on the particular subprocess that is being affected. 
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In Figure [I], we present the relative deviation of the total cross section for a par- 
ticular coupling, i.e. A' 211 , and different values of the squark masses (for convenience, 
we assume that ul and d R are degenerate). The lower cut on the dilepton invariant 
mass assumed for the calculation is 500 GeV. As far as the total cross section is 
concerned, the difference is only at the level of a few per cent, even for moderately 
large values of the coupling. However, this is somewhat misleading, since the de- 
viation from the Standard Model increases with the dilepton mass, and therefore 
the observed integrated effect strongly depends on the lower limit on the dilepton 
mass taken for the calculation. The extra contributions manifest themselves more 
markedly in the kinematic distributions, modifying them from those expected within 
the SM in three essential ways: (1) enhanced cross section for high £ + £~ invariant 
mass, (2) different lepton angular distribution in the £ + £~ rest frame and (3) different 
boost distribution of the £ + £~ system. 

For a given point in the supersymmetric parameter space, these differences are 
demonstrated in Figure [2|. The figure shows distributions in the invariant mass M 
of the dimuon pair, the difference in the rapidities of the two leptons Arj and the 
rapidity of the pair itself rj pair . Clearly, the use of the differential distributions would 
result in an enhanced sensitivity with respect to just the event rate comparison of 
Figure [l|. Note the substantial broadening of the distribution in M. In fact, rather 
than limiting ourselves to a study of a distribution in one of the three independent 
variables listed above, it is conceivable that the use of the full kinematic information, 
taking into account the correlation of the kinematic variables, would be even more 
powerful. In the next two sections we discuss how exactly we propose to do this, 
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Figure 2: The phase-space distributions for the Drell-Yan process at the LHC. We have 
demanded that the ^ have a minimum transverse momentum of 50 GeV each and be 
within the pseudorapidity range — 3 < n < 3. For the rapidity distributions, an additional 
cut (M > 500 GeV,) has been imposed. In each case, the lower curve corresponds to the SM 
and the upper curve to A^n = 0.5, niq = 800 GeV. 
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albeit with a choice of independent variables slightly different from that mentioned 
above. 



3. Maximum likelihood method 

Our goal, then, is to exploit to the full the differences in the kinematic distribution 
of the /x 4 fi~ pair caused by the presence of the ^ p interactions. Such a study is 
expected to yield a measurement (or at least constraints) in the 2-dimensional su- 
persymmetric parameter space, namely the mass-coupling plane. The problem of 
simultaneous determination of these two, viz. A' and rriq, is the classical problem of 
parameter estimation. In view of the systematic uncertainties (experimental due to 
the luminosity and theoretical due to parton distributions, ^-factors, etc.) in the 
normalisation of the signal, it is best to use the maximum likelihood method, which 
does not require a precise knowledge of the absolute size of the signal, but only that 
of its shape. 

As Figure |2| shows, the fall-off in the invariant mass distribution is indeed sub- 
stantially slower for the i? p -violating (RPV) contribution. For the other distribu- 
tions too, the difference is quite discernible. However, instead of using the variables 
of Figure ^|, we choose to work with an equivalent (independent) set, namely the 
momentum fractions x\ and x<± of the initial-state partons and the cosine of the scat- 
tering angle, cos 8. Neglecting the transverse momentum of the fi + n~ system, the 
event kinematics for a given event i is completely specified by these three variables. 
The log likelihood function is defined as: 



ln£ = ^ In F(A',m 5 , 24,4,008 0*) , (3.1) 



i=i 



with F given by: 



a(X',mg) d^ido^dcos^ 



<lCT (Y, 77^,4,4. cos 0*) . (3.2) 



The cross section cr(A', m,) in the denominator is the one obtained after imposing 
all the analysis cuts. The latter, of course, are chosen so as to maximize the RPV 
contribution in the signal. The best estimate of the true values of the two parameters 
is then given by the particular pair (X',mq) that maximizes log£ for a given data 
sample. In addition to the advantage that comes from using only the shape of 
distributions as mentioned above, this method also has the good features of not 
requiring any binning, as well as exploiting the correlations among the different 
variables optimally. We note here that the objective of avoiding the systematic 
(both theoretical and experimental) errors due to the imprecision in the knowledge 
of the absolute size of the signal, could, in principle, also be achieved by comparing 
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the size of the dimuon signal with the dielectron one. However, the latter method 
presupposes that no unknown physics effects exist in the dielectron mass spectrum 
and hence it is not completely model-independent. Furthermore, the comparison 
of two different spectra involves the compounding of errors, thereby affecting the 
accuracy adversely. 

4. Data Analysis 

For our analysis of the simulated data, we first generate events corresponding to an 
integrated luminosity of 100 fb _1 , which corresponds to a year of LHC operation 
in its high luminosity mode. The evaluation of the achievable sensitivity in the A' 
measurement requires the generation of many times the statistics expected for a 
single experiment, for a fine grid on A' for a few values of m,q. The parton-level 
generation of unit weight events without initial-state QCD showering is about ten 
times faster than the generation of full PYTH1A [26] events followed by the fast 
detector simulation for the ATLAS detector [27]. We will therefore first evaluate 
the statistical sensitivity of the experiment at the parton level. The results thus 
obtained will then be compared with the same procedure applied to fully generated 
events for a few selected points in parameter space, in order to evaluate the effect on 
the experimental sensitivity of the experimental smearing and of the introduction of 
initial-state QCD radiation. 

We select events with two isolated opposite-sign muons, satisfying the following 
two requirements: 1) P% > 10 GeV, |^| < 2.5 and 2) m M+M - > 500 GeV. The 
first criterion essentially ensures that the muons are visible in the detector. The 
invariant-mass cut, on the other hand, is motivated by the observation that the 
relative deviation of the cross section starts to become significant only at ~ 
rriq [19]. The optimal choice of the cut is determined not only by the squark mass 
under consideration, but also by the total number of events (in other words, the 
luminosity). However, rather than working with a squark mass- dependent cut, we 
choose to adopt the simpler strategy of a fixed choice for the invariant- mass cut. 
Approximately 7500 events survive these cuts for our choice of the luminosity. 

As is clear from the definition of the likelihood function, a complete knowledge 
of the kinematics of the event is necessary for its calculation. However, it must be 
borne in mind that, at a pp collider such as the LHC, it is not possible to know 
for certain the initial direction of the quark (in the qq hard scattering subprocess). 
Hence, only the absolute value of cos 6* is measurable, not the sign. Part of this 
information can, however, be recovered, by using the knowledge of x\ and x 2 and 
the fact that, in the proton, the x distribution for valence quarks is harder than for 
antiquarks. Arbitrarily labelling the proton beams "1" and "2", the difference x\— x 2 
can be inferred from the longitudinal momentum of the dimuon pair. 
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Figure 3: Difference x\ — %2 when the quark is picked up from the proton labelled "1" (left 
plot) and "2" (right plot). 

In Figure ||] we show the distribution of x\ — X2, when the quark is picked up 
from the proton labelled "1" (left plot) and "2" (right plot). It can be clearly seen 
that in most cases the quark is taken from the side with higher x. To make use of 
this information, the likelihood is built by summing the functions F calculated for 
both signs of cos#, weighted by the probability, for the given (xi, £2) combination 
that the quark is in proton "1" or "2". 

4.1 Parton level analysis 

For the sensitivity study, we generated a grid of points with different values of A' (with 
a spacing of 0.025) for four different equispaced values of the squark mass ranging 
from 300 to 900 GeV. For each point, 5 x 10 5 events with unit weight were generated. 
The likelihood was calculated in steps of A' 2 for —0.05 < A' 2 < 0.2. The inclusion of 
the unphysical range is necessary for the evaluation of the confidence interval. If the 
likelihood function has a local maximum for a positive (physical) A' 2 , we take this 
value as a measurement. Otherwise, we take the absolute maximum, even though 
it may be in the unphysical (negative) region. The cross section being a quadratic 
function of A' 2 , a secondary maximum may appear, and often does for negative A' 2 . 
In fact, even for A' 2 values large enough for the experiment to be sensitive to squark 
exchange, a small fraction of the Monte Carlo experiments can show the unphysical 
maximum to be higher than the physical one. This feature is illustrated in Figure f|, 
where we display the likelihood as a function of A 2 2 n for four experiments, assuming 
rriq = 500 GeV. The upper plots are for A 211 = 0.15, which is below the experimental 
sensitivity for this particular squark mass. In this case, even when the physical 
maximum is found, the dip between the two maxima is shallow (la ~ Aln£ = 0.5), 
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and very often the absolute maximum occurs at a negative A' 2 . The lower plots are 
for X' 2U = 0.2. The two maxima here are always well separated; the rare cases when 
the absolute maximum is in the unphysical region correspond to experiments for 
which the positive maximum is somewhat displaced from the nominal value. 
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Figure 4: Shape of the log likelihood function as a function of X^u for four different Monte 
Carlo experiments. The value of A' 211 used in the generation of events, shown as a vertical 
line, is 0.15 for the upper plots, and 0.2 for the lower plots. The assumed squark mass is 
niq = 500 GeV. The integrated luminosity is 100 fb^ 1 . 



In order to evaluate the uncertainty on the A' measurement, we generated ~ 
1500 Monte Carlo experiments for an integrated luminosity of 100 fb -1 ; for each of 
them we calculated the maximum of the likelihood function according to the above 
prescription. Since the generated statistics is only ~ 70 times the statistics for a 
single experiment, each of the Monte Carlo experiments was produced by randomly 
picking, inside the available statistics, the ~ 7500 events corresponding to one year of 
running. With this procedure each event is used for ~ 25 Monte Carlo experiments. 
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Figure 5: Distribution of the A^n value estimated through the maximization of the like- 
lihood function for a set of ^1500 Monte Carlo experiments. The upper (lower) set cor- 
responds to a squark mass of 300 (700) GeV. For each case, we show the distributions 
for three values of X' used in the generation of events: (dashed line), a value below the 
experimental sensitivity for the assumed squark mass and integrated luminosity (left plot, 
full line), and a value for which a good A measurement is achievable (right plot). The 
integrated luminosity is 100 flT 1 . 



A common behaviour is observed for all four squark masses considered. For 
values of A' 2 approximately up to the minimum in the cross section shown in Figure [l], 
the distribution of the estimated values of A' 2 is very broad, and extends up to a 
specific value of A' 2 , which depends on the mass, in such a way that the distribution 
of the measured A' 2 becomes gaussian only for A /2 values above this. This means 
that the experiments start being sensitive to the effect somewhere in this transition 
region. Figure ^| shows the results for = 300 and 700 GeV. 
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4.2 Evaluation of experimental sensitivity 

We follow here the frequentist approach of [28], where the sensitivity of an experi- 
ment is defined as: "the average upper limit that would be obtained by an ensemble 
of experiments with the expected background and no true signal". For each con- 
sidered value of the squark mass, we build, therefore, confidence belt according to 
the Neyman construction. We adopt as the auxiliary choice for the definition of the 
belt the one leading to "upper confidence intervals", defined by Eq. (2.5) in [28]. 
For an ensemble of experiments with A = 0, we then calculate the respective upper 
limits using the confidence belt thus built, and then take the average. The results are 
shown on the (m^— A') plane in Figure |6], together with the region corresponding to 
the bounds on X' 2U from low-energy processes. If no signal is present in the data, the 




200 400 600 800 1000 



m squar k (GeV) 

Figure 6: 95% sensitivity region in the rriq-X' plane for an integrated luminosity of 
100 pb^ 1 . The dashed line corresponds to the choice of the central interval in the Neyman 
construction. The analysis is performed at the parton level. The contours corresponding 
to an uncertainty on X' of 2, 5 and 10% (under the assumption of perfectly known squark 
mass) are also shown. The shaded region on the left is excluded by low energy measure- 
ments. 

result of the experiment will be the exclusion of the region in the (rrig-X') plane above 
the curve labelled '95% exclusion'. If a signal is present, the experiments will be able 
to extract a measurement of the A' couplings. The assessment of the precision of this 
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measurement depends on the available information on rriq, and on the dependence 
of the measured value of A' on the assumed squark mass. As a first, unrealistic, 
approximation we show in Figure ^| the curves in the (ra^-A') plane corresponding to 
a statistical uncertainty on A' of 2, 5 and 10% respectively, if no error is assumed on 
the squark mass. 
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Figure 7: Fractional deviation of the measured X' value from the true one as a function 
of X' for three values of the squark mass used in the likelihood calculation: the nominal 
one (dots), reduced by 3% (squares), augmented by 3% (triangles). The four plots are 
respectively for the four squark masses considered in the analysis, mg = 300, 500, 700 and 
900 GeV 



It is interesting to consider the dependence of the error estimate on the available 
prior information on the squark mass. Consider the pair-production of the said 
squark at the LHC. For the range of values of the A' couplings addressed in this paper, 
every supersymmetric event will contain two Xi decaying into two jets and a muon 
or a muon neutrino. Detailed studies have shown that, in the ATLAS experiment, it 
will be possible to reconstruct the x? peak from its decay products, and, going up the 
decay chain, to reconstruct the squark masses [29] . The statistical precision of the 
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mass measurement is essentially a function of the squark mass, and for squark masses 
below ~ 1 TeV, the error on the mass will be dominated by the systematic uncertainty 
on the reconstruction of the parton energy from the jet energy measurement. Such 
uncertainty is estimated to be approximately 3% in ATLAS [30]. Therefore, in the 
sensitivity region we have performed the A' measurement for values of the squark 
mass displaced upwards and downwards by 3% with respect to the nominal mass. 
The results are shown in Figure [7] for the four mass values considered in the analysis. 
The correlation between the A' and rriq measurement is positive, therefore a higher 
value of irtq in input yields a higher A' measurement. For the assumed value of the 
uncertainty on rriq, the additional uncertainty on A' from this effect is of order 2-3% 
for all considered mass values. One can also observe from Figure |7| that when the 
nominal mass value is used, the bias on A' from the likelihood fit is less than 1%. 



100 fb" 1 




200 400 600 800 1000 



m squark ( GeV ) 

Figure 8: la ellipses for the measurement of X' and rriq if no previous information on rriq 
is assumed for different points in the rriq-X' plane and an integrated luminosity of 100 pb^ 1 . 
The full and dashed lines correspond to the 95% exclusion region. The analysis is performed 
at the parton level. The shaded region on the left is excluded by low energy measurements. 

Finally, we consider the combined precision, in the measurement of rriq and A', 
that may be achieved if a two-dimensional likelihood is used, and no prior information 
on the squark mass is assumed. We have performed this study for a few sample points, 
and the results are shown, in Figure as la ellipses in the (rriq-X') plane, calculated 
on a sample of ~ 500 Monte Carlo experiments. The resolution in mass is of the 
order of few tens of GeV, and quickly degrades when A' approaches the edge of the 
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m g (GeV) 


A' 


a(X')(l) 


a(X')(2) 


a(m d ) (GeV) 


p 


300 


0.175 


0.005 


0.014 


43 


0.913 


500 


0.200 


0.010 


0.027 


114 


0.947 


500 


0.250 


0.007 


0.020 


63 


0.948 


500 


0.275 


0.006 


0.017 


46 


0.951 


500 


0.300 


0.005 


0.014 


38 


0.949 


700 


0.250 


0.013 


0.036 


148 


0.859 


700 


0.300 


0.009 


0.028 


98 


0.951 


700 


0.350 


0.007 


0.020 


56 


0.954 


700 


0.400 


0.005 


0.016 


40 


0.964 



Table 1: Statistical errors on the X' measurement for a few sample points. The value in the 
third column [a(X')(l)J is the uncertainty assuming that the squark mass is known with zero 
error. The last three columns give the parameters of the la ellipses for the two-dimensional 
likelihood on (rriq-X'), where p is the correlation coefficient for uig and A'. 



sensitivity region. Given the correlation between the A' and m q measurements, the A' 
measurement is degraded accordingly. This statement is better quantified in Table 1, 
where we give the parameters of the ellipses corresponding to log(£ max ) — 1/2 for the 
studied points, compared with the resolution on A' if the squark mass is known. In 
particular, the parameter p, which measures the correlation of the two variables, is 
~ 1, meaning that the variables are fully correlated. Thus, a combined measurement 
of squark mass and A' is possible with this analysis for a significant fraction of the 
parameter space. 



5. Analysis of fully generated events 

The analysis in the previous section was performed assuming that the momentum of 
the fi + n~ system has no component transverse to the beam axis, and that the muon 
momenta are perfectly measured in the detector. 

In order to perform a more realistic evaluation of the precision of A' measurement, 
achievable in a real experiment, the matrix elements for the squark exchange process 
have been inserted into the PYTHIA [26] event generator as an external process, 
and full events have been generated, including the full PYTHIA machinery for QCD 
showering from the initial-state quarks, and for the hadronization. The events thus 
generated have been passed through the fast simulation of the ATLAS detector [27], 
including, in particular, a very detailed parametrization of the resolution of the muon 

momentum, measurement. , . r , r , , , r r , , 

lhe values of xi, x 2 have been interred from the measured tour-momenta ot the 
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Figure 9: Left plot: correlation of the reconstructed value of log(xi) with the generated 
value for a sample of events passing the analysis requirements. Right plot: distribution of 
the ratio of the reconstructed and generated value of x\ . The distributions are given for a 
sample of events passing the analysis cuts. 



detected muons, according to the formulae: 



2 Pf»~ 



Xi 



x 2 , m u + 



MA* 



X1X2S 



For the evaluation of cos# we use the Collins-Soper convention [31], which consists of 
an equal sharing of the system transverse momentum between the two quarks. 

We show, in Figure 0, the two-dimensional correlation between the generated and the 
reconstructed values for the variable x\ as well as their ratio. Only events that pass 
the analysis criteria enter these plots. For x\, the RMS deviation is ~ 10%, and is 
completely determined by the muon momentum resolution. A comparable resolution 
is obtained for | cos6*| (Figure [TUp, but with significant tails due to the presence of a 
non-zero transverse momentum of the fi + fi~ system. 





0.25 0.5 0.75 

|cos9|(true) 




0.5 1 1.5 

|cos9|(measured)/|cos9|(true) 



Figure 10: As in Figure but for \ cos I 
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Samples of events were generated for four squark masses, 300, 500, 700 and 
900 GeV, and for three values of A' for each mass, chosen in each case in an interval 
of 0.1 in A' 2 , approximately covering the transition region where this analysis be- 
comes sensitive to the RPV contribution for an integrated luminosity of 100 fb _1 , as 
described in the previous section. We then applied to these samples the unbinned 
likelihood analysis, as for the parton-level events. The evolution of the likelihood 
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Figure 11: Distributions of measured X' 2 for a set of Monte Carlo experiments for different 
values of X' used for the generation of events. The full-line histograms are the parton-level 
distributions while the dashed histograms correspond to fully generated events. 

distributions with A' is nicely reproduced. In particular, the threshold behaviour of 
the sensitivity is still present. The distribution of the measured A' is indeed domi- 
nated by the statistical performance of the maximum likelihood estimator, and not 



by the experimental resolution. The results of the analysis are shown in Figure [I I 
for uiq = 300 GeV. Therefore, the sensitivity limit calculated at parton level, shown 
in Figure |], applies for the full event simulation as well. A difference is observed only 
for A' values for which the distribution of the measured A' for a set of Monte Carlo 
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experiments is approximately gaussian. In this case, the presence of the experimental 
smearing induces a moderate increase in the A' measurement resolution. Figure [12| 
shows a comparison of the achievable parton-level resolution with the corresponding 
experimental resolution. The four continuous lines give the estimated parton-level 
resolution for the four considered masses. Superimposed as data points are the reso- 
lutions obtained for the experimental simulation. A difference is observed mainly for 
the lower A' values, where the experimental simulation exhibits higher non-gaussian 
tails than the parton-level simulation, thus yielding a 10-20% higher RMS deviation. 
For higher A' values, the resolutions are essentially identical. We can thus conclude 
that the curves shown in Figure || are only marginally affected by the experimental 
uncertainties. This result is obtained under quite pessimistic assumptions. In fact, 
we use as a definition of F and its integral in Eq. (|3.2Q the pure leading order formulas 
in Eq. ( |2.2j ), without any attempt at parametrizing the fi + fi~ transverse momentum 
distribution. 




0.1 0.2 0.3 0.4 

X 

Figure 12: Curves of a(X') / X' as a function of X' for an integrated luminosity of 100 fb^ 1 
and for four different squark masses. The full lines give the resolution obtained at parton 
level, and the points are the results of the full event simulation. 
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In the analysis of real data, we also need to consider the backgrounds from non- 
Drell-Yan processes. The dominant background processes, for the invariant mass 
cut applied in this analysis, will be bb, it, and WW production. For bb production 
the muons are not isolated, and we can apply lepton isolation, which consists of 
requiring an energy deposition of less than 10 GeV, not associated with the lepton 
in a pseudorapidity-azimuth (77— 0) cone of opening AR = 0.2 around the lepton 
direction. 

Without any additional cuts, for an integrated luminosity of 100 fb _1 , the back- 
grounds from it and WW amount to ~ 1400 and ~ 500 events respectively, for a 
signal of ~ 7400 events. The it background can be strongly reduced by vetoing 
6-jets. We assume a tagging efficiency of 50% for a rejection factor of 100 (10) on 
light quark (charm) jets respectively. This is appropriate for the high luminosity 
run. We also veto jets tagged as 6-jets with Pt > 30 GeV. The it background is 
thus reduced to ~ 600 events, with negligible effect on the signal. Both for the it 
and WW backgrounds, the events will have real E™ lss from escaping neutrinos, and 
a further reduction can be achieved by vetoing on high E™ tss events. This require- 
ment, though, has a significant effect on the kinematics of signal events. In fact we 
are considering here high energy muons, for which the error in momentum measure- 
ment induces an instrumental imbalance in the vector sum of the lepton momenta, 
which grows with increasing momentum. Therefore the acceptance of any kinematic 
cut applied must be convoluted in the test function used for the likelihood in order 
to obtain the correct result. A complementary approach consists in accepting the 
relatively high level of background, and incorporating the background shape into the 
likelihood function. This should be possible with high precision, since the considered 
backgrounds yield twice as many e/i events as /x/x events, allowing an estimate of 
the background level from the data themselves. The two approaches have different 
systematic uncertainties, and can be used in parallel, thus providing a double check 
on the result. 

6. Systematic uncertainties 

In the previous section, we have studied in detail the main sources of experimental 
uncertainty, coming from the imperfect measurement of the event variables. Other 
possible sources of experimental error are the uncertainties in the muon energy scale, 
the linearity for high-energy muons and the acceptance evaluation. The assessment of 
the effect of these uncertainties would require a detailed detector simulation, which is 
outside the scope of this work. We should however remark that the analysis described 
in the previous section displays little sensitivity to the details of the modelling of the 
experimental resolution and acceptance. In fact, in building the likelihood, the cru- 
cial factor is the normalization integral in the denominator of Eq.( |3.2| ). That integral 
is calculated with a sharp cut on the generated lepton pseudorapidity and invariant 
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mass for events with no P^t p . The data selection for the particle level analysis is 
instead applied on leptons that have been smeared according to the detector reso- 
lution. Therefore, the acceptance of the cuts, especially the one on m^+^~ edges, is 
reproduced in a very approximate way. Notwithstanding this fact, we observe a good 
agreement between the parton-level analysis and the particle level one. This gives 
us confidence that the additional experimental uncertainty on the A' measurement 
resolution would, at worst, only marginally affect the results of our analysis. More 
important would be the effect of a non-linearity in the lepton energy measurement, 
which could simulate a deviation from the invariant-mass distribution predicted for 
the Standard Model. A control on the linearity at the per cent level will be needed 
for this analysis. 

We also need to consider the theoretical uncertainties in the likelihood calcula- 
tion. The likelihood function is built by weighting real events according to a theoreti- 
cal cross section formula. Any discrepancy between the theoretical formula employed 
and reality will induce an uncertainty on the measurement of A'. Two main sources 
of uncertainty can be identified: 1) the presence of higher-order correction to the 
processes and 2) the parton distribution function (PDF) for the proton. We consider 
here the most important QCD higher-order corrections coming from radiation from 
initial-state quarks, which impart a transverse momentum to the lepton-lepton sys- 
tem. From the above discussion of the analysis, it is clear that the results show very 
little sensitivity to the detailed modelling of the transverse momentum distribution 
of the dilepton system. In fact, even for the fully generated events the likelihood was 
built from the pure leading order formulas, whereas the events are generated with 
the full PYTHIA machinery for initial-state radiation. Therefore, the experimental 
error quoted in the previous section implicitly includes very pessimistic assumptions 
on our ability to model this effect. In a real experiment, a more realistic theoretical 
modelling will probably be used to build the likelihood. 

The second source of uncertainty is more important, since it will affect the shapes 
of the distribution. All the events were generated with the CTEQ4L PDF set. We 
have evaluated the effect of an uncertainty on PDF parametrization by performing 
the likelihood using sets of PDFs different from the ones used for event generation. 
We chose the PDFs labelled as CTEQ4A1 and CTEQ4A5, which span an extreme 
range between 213 and 399 MeV for the QCD parameter A' 4 l The value of A^ for 



the set CTEQ4L is 236 MeV. We show, in Figure [13], the fractional deviation of the 
measured A' from the true one for the default structure functions CTEQ4L as well 
as for two other possible choices. The deviation for CTEQA1 is small, as expected 
from the small difference in AM. For the set CTEQA5 the fractional difference 
is approximately 10% for the lowest values of A' considered, and it decreases with 
increasing A' to a few per cent. We can tentatively conclude that the uncertainty on 
the measurement coming from the choice of PDFs will be at the few per cent level, 
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Figure 13: Fractional deviation of the measured X' value from the true one as a function 
of X' for three different choices of structure functions in the calculation of the likelihood. 
The events were generated with CTEQ4L. The upper plot is for niq = 500 GeV, the lower 
one for niq = 700 Ge V 

comparable to the statistical one. 



7. Effect of the $ p couplings on the decays of the squarks 

The above discussion establishes that the dimuon production at the LHC has enough 
sensitivity to $ p couplings to admit their measurement with a fair accuracy. Our 
analysis shows that a Ijtp coupling of the order of the gauge coupling can be measured 
quite accurately even for rather large squark masses, up to 700-900 GeV. Clearly, 
squarks in this mass range have a substantial production cross section at the LHC. 
So this naturally brings us to the question of whether one could detect or even 
measure such couplings in the decays of these squarks. A discussion of this issue is 
clearer if we review a few facts about the possibilities of measuring the ty p couplings 
in collider experiments in general. These are, of course, strongly dependent on the 
value of the Ifi p coupling being probed. In the worst-case scenario, it may be so 
small that the only effect it has is to cause the decay of the LSP and that too 
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outside the detector, so that it is mistaken for a stable non-interacting particle. 
In this situation, we will altogether miss the phenomenon of R-parity violation at 
the colliders. For intermediate values of these couplings, one will get to a region 
where a significant fraction of the LSPs will decay inside the inner cavity of the 
LHC detector and in that case one might be able to measure the strength of these 
couplings by detecting the displaced vertices inside the detector. To the best of our 
knowledge, a detailed experimental study to map the parameter space in this case 
is not currently available. Then comes the range of couplings that are large enough 
to cause a prompt decay of the LSP, giving rise to striking final states, but are still 
not large enough to cause single-sparticle production or significantly affect decays 
of sparticles other than the LSP. In this range of values the distinctive final states 
caused by the decay of the LSP due to the fy p couplings can provide evidence for the 
existence of Ijt p , but cannot give any information whatsoever on their size. And finally 
comes the region of the even larger Ifi p couplings that we have considered. Here, one 
is sensitive to the effects of these couplings via virtual exchanges of sparticles on 
scattering processes as well as decays of sparticles, other than the LSP, caused by 
them. We have demonstrated in the above work the feasibility of 'measuring' such 
a Ijt p coupling through the contribution of virtual squark exchanges to the dimuon 
production. This offers perhaps the cleanest way of measuring the ]fi p couplings with 
an adequate control of the experimental and theoretical systematic uncertainties, as 
the only input required for this is the squark mass. As already mentioned, for values 
of the squark masses under consideration, their production cross sections at the LHC 
are substantial. Hence, the effect of the rf, p couplings on the decays of the squark 
produced via the strong interactions needs to be studied. A particularly attractive 
case will be when the Ifi p decay of the squark competes with an /^-conserving decay 
with an identifiable signature. A comparison of the two can then give a direct 
measurement of the ]fi p coupling. A promising example for the couplings we have 
considered is the ]fi p decay t — > fj,+ jet. To study this issue, we will need to consider 
both the total production cross sections and the branching fractions into the various 
channels. 

Squark pair-production has been analysed at length in various papers and de- 
pends crucially on the gluino mass. A crude summary of these results relevant to 
our discussion is that "for any squark lighter than approximately 800 GeV, the pair- 
production cross section at the LHC is large enough to lead to discovery". For the 
most part, the reach in mass and possibility of discovery will be independent of 
whether .R-parity is violated or not. The measurement of the relevant coupling, or 
even the establishment of its identity, would, however, require that at least two chan- 
nels be measurable with a relatively high degree of accuracy. One of these should 
be i?-parity-violating, while the other should be /^-conserving. For top-squarks, the 
presence of leptons and 6-jets in the final state provides a handle for the identifica- 
tion of the exclusive decays of interest. The relevant .R-conserving decay channels are 
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Figure 14: Contours for the branching fraction of a 300 GeV stop (t^) into the R-parity- 
violating channel. The unification relation between the gaugino masses has been assumed. 
The value of the R-violating coupling (X^i = 0.16) is the minimum necessary for a deter- 
mination X' to an accuracy of 5% (see Figure ^). The left and right panels correspond to 
two different values o/tan/3. 

q — ► qg, qx®, q'xf, q'W. Each of these channels could lead to a cascade process culmi- 
nating in a Ifip decay. We are, however, concerned with the Ifi p decay of the original 



squark. In Figure 14, we display ^-branching fraction contours (into /i + + d) for the 
decay of a 300 GeV stop. The assumed value of A' 231 is the least that would permit 
a detection with a statistical accuracy of 5%, if any dilution in the efficiency of the 
analysis from surviving backgrounds is disregarded. We see that barring very light 
charginos, the ]fi p mode tends to dominate. An overwhelmingly large IjL p branch- 
ing fraction, while leading to spectacular signals, is hardly amenable to a precision 
measurement of the coupling. Similar arguments hold for a very small F/L p branching 
fraction. The best hope is therefore for moderate values of M 2 and fi. It is thus clear 
that, at the level of branching ratios, stop decays would provide a sufficient number 
of events for the measurement of the I/L p coupling over a significant parameter space. 
A detailed experimental analysis, outside the scope of this work, is however needed 
to ascertain whether it will be possible to isolate the final states produced by the two 
relevant processes, with adequate efficiency, and with a level of purity that would 
permit a reasonably accurate determination of the relevant branching fractions. 

8. Conclusions 

To summarize, we demonstrate that an analysis of the Drell-Yan process at the LHC 
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detectors would be a significant tool in the task of probing the parameter space in an 
i?p-violating supersymmetric model. While the deviation in the total cross section 
is at best a few per cent, even for large values of the relevant R p - violating coupling 
constant, the differential distributions are much more discriminating, the distribution 
in the dilepton invariant mass proving particularly useful. By adopting the maximum 
likelihood method we could maximize the sensitivity of the measurement, avoiding 
at the same time the normalization uncertainties due to structure functions as well 
as higher-order corrections, etc. 

Working, for definiteness, with dimuon production in the ATLAS detector, we 
showed that, for a wide range of the parameter space, i? p -violating supersymmetry 
would be amenable to discovery through this process. Even more importantly, a 
measurement of the squark mass as well as the Tji v coupling would be possible with 
relatively low errors: 2-3% for the coupling and somewhat larger for the mass. These 
errors could be reduced even further if additional information about the squark mass 
were to become available from other measurements. We found the systematic errors 
to be small. The analysis proved to be quite robust and the results gleaned from 
fully generated events were very similar to those obtained from a parton-level study. 
Since for the range of squark masses to which our investigation is sensitive, they can 
be pair-produced copiously at the LHC, we also further looked into the possibility 
of determining the same F/L p coupling, for which we analysed the sensitvity using the 
Drell-Yan process, through the decays of the squarks. 

Finally, let us add that even though we have confined ourselves to dimuon pro- 
duction in a particular theory (namely R p - violating supersymmetry), a similar anal- 
ysis could be carried out for dielectron production equally well. In general such an 
analysis can be used to study effects of any alternate theories of physics, beyond the 
SM, which affect the Drell-Yan process. Our analysis demonstrates that such studies 
would complement very well the more direct methods for new particle search. 
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